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(N 1 Abstract 
>V 

We consider control systems of the type x = Ax + a(t)bu, where u € R, (A, b) is a controllable 
pair and a is an unknown time- varying signal with values in [0, 1] satisfying a persistent excitation 
condition i.e., f^ +T a(s)ds > \i for every t > 0, with < fi < T independent on t. We prove 
that such a system is stabilizable with a linear feedback depending only on the pair (T, fi) if the 
eigenvalues of A have non-positive real part. We also show that stabilizability does not hold for 
arbitrary matrices A. Moreover, the question of whether the system can be stabilized or not with 
, an arbitrarily large rate of convergence gives rise to a bifurcation phenomenon in dependence of 

C$ . the parameter fi/T. 



m ■ 1 Introduction 

> 

CN ■ The present paper is a continuation of [9] , where the study of control linear systems subject to scalar 

' persistently excited PE-signals was initiated. The general form of such systems is given by 
O ; 

x = Ax + a{t)Bu , (1) 

o : 

J> , where x G R", u £ R m and the function a is a scalar PE-signal, i.e., a takes values in [0, 1] and 
^ there exist two positive constants fj,, T such that, for every t > 0, 

rt+T 

a(s)ds > fi. (2) 



*The work was in part carried out while the first author was working as Marie Curie Fellow at the Department 
of Mathematics and Statistics, University of Kuopio, Finland, supported by the European Commission 6th frame- 
work program "Transfer of Knowledge" through the project "Parametrization in the Control of Dynamic Systems" 
(PARAMCOSYS, MTKD-CT-2004-509223). The second author was partially supported by the grant ArHyCo of the 

Agence Nationale de la Recherche. 

' Laboratoire des Signaux et Systemes, Supelec, 3, Rue Joliot Curie, 91192 Gif s/Yvette, France and Universite 

Paris Sud, Orsay, chitour@lss.supelec.fr 

* Institut Elie Cartan, UMR 7502 INRIA/Nancy-Universite/CNRS, BP 239, Vandoeuvre-les-Nancy 54506, France 
and CORIDA, INRIA Nancy - Grand Est, France, mario.sigalotti@inria.fr 



1 



Given two positive real numbers \i and T with fi < T, we use Q(T,fi) to denote the class of all PE 
signals verifying (2). 

In (1), the PE-signal a can be seen as an input perturbation modeling the fact that the instants 
where the control u acts on the system are not exactly known. If a only takes the values and 1, then 
(1) actually switches between the uncontrolled system x = Ax and the controlled one x = Ax + Bu. 
In that context, the PE condition (2) is designed to guarantee some action on the system. (For 
a more detailed discussion on the interpretation of persistently excited systems and on the related 
literature, see [9].) 

Our main concern will be the global asymptotic stabilization of system (1) with a constant linear 
feedback u = —Kx where the gain matrix K is required to be the same for all signals in the 
considered class G(T,fi) i.e., K depends only on A,b,T,fj, and not on a specific element of G(T,fi). 
We refer to such a gain matrix K as a (T, /i)-stabilizer. It is clear that (^4, B) must be stabilizable 
for hoping that a (T, /i)-stabilizer exists and we will suppose that throughout the paper. Moreover, 
the stabilizability analysis can be reduced to the controllability subspace and thus to the case where 
(^4, B) is controllable. 

The questions studied in this paper find their origin in a problem stemming from identification 
and adaptive control (cf. [3]). Such a problem deals with the linear system x = —P(t)u, where the 
matrix P(-) is symmetric non-negative definite and plays the role of a. If P = /, then u* = x trivially 
stabilizes the system exponentially. But what if P(t) is only semi-positive definite for all t? Under 
which conditions on P does u* = x still stabilize the system? The answer for this particular case, 
can be found in the seminal paper [13] which asserts that, if x G R n and P > is bounded and has 
bounded derivative, it is necessary and sufficient, for the global exponential stability of x = —P(t)x, 
that P is also persistently exciting, i.e., that there exist /i,T > such that 

j^ +T ?P(s)Zds>n, (3) 

for all unitary vectors £ £ R n and all t > 0. Therefore, as regards the stabilization of (1), the notion 
of persistent excitation seems to be a reasonable additional assumption on the signals a. 

Let us recall the main results of [9]. We first addressed the issue of controllability of (1), uniformly 
with respect to a G G(T, /i). We proved that, if the pair (A, B) is controllable, then (1) is (completely) 
controllable in time t if and only if t > T — \x. We next focused on the existence of (T, ^-stabilizers. 
We first treated the case where A is neutrally stable and we showed that in this case the gain K = B T 
is a (T, /i)-stabilizer for system (1) (see also [3]). Note that in the neutrally stable case K does not 
depend on T and /i. We next turned to the case where A is not stable. In such a situation, even 
in the one-dimensional case, a stabilizer K cannot be chosen independently of T and /U. In [9], we 
considered the first nontrivial unstable case, namely the double integrator x = J2X + abou, where 
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J2 denotes the 2x2 Jordan block, the control is scalar and 60 = (0, 1) T ■ We showed that, for every 
pair (T, //), there exists a (T, /^-stabilizer for x = J2X + ab^u, a € Q{T,p). 
In this paper, we restrict ourselves to the single-input case 

x = Ax + a(t)bu, u £ R, a£Q(T,p), (4) 

and we provide two sets of results. The first one concerns the stabilizability of (4). Given two 
arbitrary constants p and T with < p <T, we prove the existence of a (T, /i)-stabilizer for (4) when 
the eigenvalues of A have non-positive real part. The second set of results concerns the possibility 
of obtaining an arbitrary rate of convergence once stabilization is achieved. We essentially focus on 
the two-dimensional case and we point out an interesting phenomenon: there exists S (0, 1) so 
that, for every controllable two-dimensional pair (A, b), every T > and every p £ (0, p*T), the 
maximal rate of convergence of (4) is finite. Here maximality is evaluated with respect to all possible 
(T, //)-stabilizers. As a consequence, we prove the existence of matrices A (e.g., J2 + AW2 with A 
large enough) such that for every T > and every p S (0, p*T), the PE system (4) does not admit 
(T, //)-stabilizers. The latter result is rather surprising when one compares it with the following two 
facts: let p G (0, 1]; (i) given a sequence (a n ) nS N with a n £ Q(T n ,pT n ) and lim n ^ +00 T n = 0, all its 
weak-* limit points a* take values in [p, 1] (see Lemma 2.5) and (ii) the two-dimensional switched 
system x = J2X + a+bou can be stabilized, uniformly with respect to a+ £ L°°(R>o, [p, 1]), with an 
arbitrary rate of convergence. The weak-* convergence considered in (i) is the natural one in this 
context since it renders the input-output mapping continuous. 

Let us briefly comment on the technics used in this paper. First of all, it is clear that the notion of 
common Lyapunov function, rather powerful in the realm of switched systems, cannot be of (direct) 
help here since, at the differential level, one can evolve with an unstable dynamics x = Ax, when 
a = takes the value zero. More refined tools as multiple and non-monotone Lyapunov functions 
(see, e.g., [1, 2, 7, 10, 14, 16]) do not seem well-adapted to persistently excited systems, at least for 
what concerns the proof of their stability. It seems to us that one must rather perform a trajectory 
analysis, on a time interval of length at least equal to T, in order to achieve any information which 
is uniform with respect to a S Q(T,p). This viewpoint is more similar to the geometric approach 
to switched systems behind the results in [4, 5, 6]. As a second consideration, notice that point 
(i) described above, which is systematically used in the paper, presents formal similarities with 
the technique of averaging but is rather different from it, since no periodicity nor constant-average 
assumption is made here. Moreover, for a given persistently excited system, T is fixed and thus it 
does not tend to zero. 

The paper is organized as follows. In Section 2 we introduce the notations of the paper, the basic 
definitions and some useful technical lemmas. We gather in Section 3 the stabilizability results for 
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matrices whose spectrum has non-positive real part. Finally, the analysis of the maximal rates of 
convergence and divergence is the object of Section 4. Since many of our results give rise to further 
challenging questions, we propose in Section 5 several conjectures and open problems. 

2 Notations and definitions 

Let N denote the set of positive integers. Given n and m belonging to N, we use nxm to denote 
the n x m matrix made of zeroes, M n (R) the set of real-valued n x n matrices, and Id n the n x n 
identity matrix. We also write n for nX i, cr(A) for the spectrum of a matrix A G M n (R), and 9ft (A) 
(respectively, 3f(A)) for the real (respectively, imaginary) part of a complex number A. 

Definition 2.1 (PE signal and (T, /i)-signal) Let ji and T be positive constants with n < T. A 
(T, /i)-signal is a measurable function a : R>o — ► [0,1] satisfying 

t+T 

a(s)ds > ll , Vi G R>o • (5) 

We use Q(T,/j.) to denote the set of all (T, li)- signals. A PE signal is a measurable function a : 
R>o — > [0, 1] such that there exist T,/j, for which a is a (T, fi) -signal. 

Definition 2.2 (PE system) Given two positive constants fx and T with fx <T and a controllable 
pair (A,b) G M n (R) x R n , we define the PE system associated to T,li,A, and b as the family of 
linear control systems given by 

x = Ax + aub, a G G(T, ll). (6) 

Given a PE system (6), we address the following problem. We want to stabilize (6) uniformly with 
respect to every (T, /i)-signal a, i.e., we want to find a vector K G R n which makes the origin of 

x = (A - a(t)bK T )x (7) 

globally asymptotically stable, with K depending only on A, b, T and ll. 

More precisely, referring to x{- ;to,XQ,K, a) as the solution of (7) with initial condition x(to> ^0; ^0) K, a) = 
xo, we introduce the following definition. 

Definition 2.3 ((T, /i)-stabilizer) Let n and T be positive constants with Li < T. The gain K is 
said to be a (T, /i)-stabilizer for (6) if (7) is globally asymptotically stable, uniformly with every (T, a)- 
signal a. Since (7) is linear in x, this is equivalent to say that (7) is exponentially stable, uniformly 
with respect to a G Q(T,ll), i.e., there exist C, 7 > such that every solution x(- ;to,XQ,K,a) of (7) 
satisfies 

||x(t;to,x ,iiC,a)|| < Ce~ (t_to)7 ||a;o||, Vt > to- 
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The next two lemmas collect some properties of PE signals. 
Lemma 2.4 1. If a(-) is a (T, p) -signal, then, for every t$ > 0, the same is true for a(to + •). 

2. IfO<p'<pandT>0 then g(T, pT) C G(T, p'T). 

3. For 7] G (0,fj), G(T,p) C Q(T + rj, p) n G(T - rj, p - rj). 
4- IfT>r>0 and p>0, then G{t, P t) C G(T,(p/2)T). 

5. For every < p' < p there exists M > such that for every T > Mr > one /ias G(t, pr) C 

a(r, P 'r). 

Proof. We only provide an argument for points 4 and 5. Fix t > 0, T > t, p > and a G C/(t, pr). 
Let Z be the integer part of T/r. Since Z > max(l, T/t — 1), then f* +T a(s)ds > Zpr > max(r, T — 
r)p > Tp/2. For pf G (0,p) and T/t large enough, one has max(r, T — r) > (p' / ' p)T and so 
Jl +T a{s)ds> p'T. M 

Let 
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Recall that an element / of L°°(R>o, [0, 1]) is the weak-* limit of a sequence (fk)keN of elements 
of L°°(R> , [0, 1]) if, for every g G ^(R^R), 



/•oo 

f(s)g(s)ds = lim / f k (s)g(s)ds. 
k ^°° Jo 



(8) 



It is well known that L°°(R>o, [0, 1]) endowed with the weak-* topology is compact (see, for instance, 
[8]). Hence, each G(T,p) is weak-* compact. Unless specified, limit points of sequences of PE signals 
are to be understood as limits of subsequences with respect to the weak-* topology of -L°°(R>o, [0, 1]). 

Lemma 2.5 Let (c/ n )) n eN md (f n )neN be, respectively, a sequence of (T, p)- signals and an increas- 
ing sequence of positive real numbers such that lirrin^oo v n = oo . 

1. Define a n as the (T/u n , p/u n )- signal given by a n (t) = a^ n \v n t) fort > 0. If a* is a limit point 
of the sequence (a n )neN> then a* takes values in [p/T,l] almost everywhere. 

2. Let jo G {0, 1} and h G N. Let ujj, j = jo, . . . , h, be real numbers with u>j = if and only if 
j = and {iwj} ^ for j ^ I. For every t >0, let 

/ 1 ^ / , \ 

/ e UlAo % \ 

if jo = or v(t) = : if j = 1. 



v(t) 



e ulAot bo 



/ 



V 



a ui h Aot 



bo j 



For every signal a and every t > 0, define 

a c (i) = a(t)v(t)v(t) T . (9) 



Then a is a time- dependent non-negative symmetric (2h + 1 — jo) x (2h + 1 — jo) matrix with 
a c < ld2h+l-j and there exists £ > only depending on T, fi and ujj , . . . ,u>h such that, for 



every t > 0, 

J Cf {T)dT > £ld 2h+l _ 30 . (10) 

Define, moreover, ct^(t) = (a^) c '{y n t) for every t > and every n G N. If is a limit 
point of the sequence (a„) nS N for the weak-* topology of L ao (R>o, M2/ l +i_ JO (R)), then oq > 
(£/T)ld2h+i-j almost everywhere. 

Proof. Let us first prove point 1. Let a* be the weak-* limit of some sequence (ot nk )k>i- For every 
interval J contained in R>o of finite length |J| > 0, apply (8) by taking as g the characteristic 
function of J. Since each a nk is a (T/i/ nh ,ii/u nk )-signal, it follows that 

\JWn k \ M 



-— [ a*(s)ds= lim — f a n , k (s)ds > liminf — p — 2 

\J\Jj fc-»oo \ J\ Jj fc->oo \J\V nk 



T 



where !(■) denotes the integer part. Since a* is measurable and bounded, almost every t > is a 
Lebesgue point for a*, i.e., the limit 

1 !' t+£ 
lim — / a.As)ds 
£ -o+ 2e J t _ e ' 

exists and is equal to a*(£) (see, for instance, [15]). We conclude that, as claimed, a* > fi/T almost 
everywhere. 

For the first part of point 2 fix t > and notice that the map 

r-t+T 

C I 



a i— ► J a {s)ds 

is continuous with respect to the weak-* topology and takes values in the set of non-negative sym- 
metric matrices. 

We claim that all such matrices are positive definite. Assume by contradiction that there exist 
a G G(T,n) and x G K 2h+1 ~^ \ {0 2h+1 _ jo } such that f* + XQa c (s)xods = 0. Then, for almost 
every s G [t, t + T], we would have a(s)xQv(s) = 0. Since a(s) ^ for s in a set of positive 
measure, we deduce that the real-analytic function Xqv(-) is identically equal to zero. Let A% = 
diag(l,wii4o, . . .,uj h A ) if jo = or = diag(wiA , • • -,u h A ) if jo = 1. Then x%(A$yv(0) = 
for every non-negative integer j. The contradiction is reached, since (Aq,v(0)) is a controllable pair 
and x / 2 h+i-j - 
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Then, by weak-* compactness of G(T, fi), we deduce the existence of £ > independent of a such 
that (10) holds true. The independence of £ with respect to t follows from the shift-invariance of 
Q(T,ji) pointed out in Lemma 2.4. 

The second part of point 2 follows from the same argument used to prove point 1, noticing that, 
for every t > 0, 



a£ (r)dT > — Id 2 / l +i- 

t v n 



JO- 



3 Spectra with non-positive real part 

We consider below the problem of whether a controllable pair (^4, b) gives rise to a PE system that 
can be (T, /^-stabilized for every choice of ji and T. We will see in Section 4 that this cannot be 
done in general. The scope of this section is to study the case in which each eigenvalue of A has 
non-positive real part. 

The first step is to consider the special case of the n- integrator. Let J n £ M n (R) be defined as 

/o 1 o\ 

1 ••• 



J ri 







1 



Theorem 3.1 Let A = J n and b = (0, . . . , 0, 1) T E R n . Then, for every T, fi with T > fi > there 
exists a (T, fj,)- stabilizer for (6). 



Proof. In the special case of the n-integrator system (7) becomes 



xj = Xj+i, for j = 1, . . . , n - 1, 
x n = -a(t)(kixi H h k n x n ) , 



(11) 



where K = (fci, . . . , k n ) . 
For every v > 0, define 

As done in [9] in the case n = 2, one easily checks that, in accordance with 

vD n i/JnT)ri,v — Jni D n u b — 6, 



(12) 



(13) 
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the time-space transformation 



z„(i) = D-] v x(yt) , Vt>^, (14) 



of the trajectory x(-) = x{- ; to, xq, K, a) satisfies 



— x u (t) = J n x u (t) - a(vt)vbK D njU x u (t), 



that is, 



x u(-) = x(-;to/v,D n *xo,vD n>u K, a(v ■)). (15) 



As a consequence, (11) admits a (T, /Li)-stabilizer if and only if it admits a (T/v, /x/^)-stabilizer. More 
precisely, K is a (T, /i)-stabilizer if and only if vD n ^ v K is a (T/i/, /x/z^)-stabilizer. 
Let us introduce, for every gain K, the switched system 

Xi = , for j = 1, . . . , n — 1, 

J J+ ^G J L 00 (R>o,[///T,l]). (16) 

X n = —a.±(t)(klX\ H h/Cn^n), 

Recall that (16) is said to be globally uniformly exponentially stable as a switched system if the 
origin is globally exponentially stable, uniformly with respect to a* G L°°(R>o, \p/T, 1]), for the 
dynamics of (16). (For this and other notions of stability of switched systems see, for instance, [12].) 
For every K such that k\ / 0, define X\ = k\X\ + • • • + k n x n , X2 = k\X2 + • • • + k n -\x n , . . . , 
X n = k\x n . The global uniform exponential stability of (16) is clearly equivalent to that of 

Xj = X j+1 - aJcjXi, j = l,...,n, a*(t) G \p/T, 1], (17) 

where kj = k n+ \^j and, by convention, X n+ \ = n . 

It has been proven in Gauthier and Kupka [11, Lemma 4.0] (where the result is attributed to 
W.P. Dayawansa), that there exist K G R n , a scalar 7 > and a symmetric positive definite n x n 
matrix S such that 

(J n - aK(l, 0, . . . , 0)) T S + S(J n - aK{l, 0, . . . , 0)) < -jld n , (18) 

for every (constant) a G \p/T, 1]. 

Hence, there exist a gain K G R" such that (16) is globally uniformly exponentially stable and 
a positive definite matrix 5' such that the quadratic Lyapunov function V(x) = x T S'x decreases 
uniformly on all trajectories of (16). In particular, there exists a time r such that every trajectory 
of (16) starting in B\ = {x G R n | V(x) < 2} at time lies in B\ = {x G R n | V(x) < 1} for every 
time larger than r. 

We claim that, for some v > 0, every trajectory of x = (A — a u (t)bK T )x with initial condition 
in B\ and corresponding to a (T/u, fi/u) -signal a u stays in B\ for every time larger than It. (In 
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particular, by homogeneity, K is a (T/u, /i/^)-stabilizer and thus v~^D~\K is a (T, /i)-stabilizer.) 
Assume, by contradiction, that for every I G N there exist xo,z £ -S^) ^ £ [2t, 4r] and a; G Q(T/l, fi/l) 
such that 

x(i z ;0,x ,J, A',a ? ) for every I G N. (19) 

By compactness of B\ x [2r, 4r] and by weak-* compactness of L°°(R>o, [0, 1]), we can assume that, 
up to extracting a subsequence, xqi — > xo,* G S^, — ► G [2r, 4r] and a; converges weakly-* to 
a* G L°°(R>o, [0, 1]) as I goes to infinity. Then x(ti;0, xoj, K, ai) converges, as I goes to infinity, 
to x(t+; 0, x*, X, a*) (see [9, Appendix] for details). Since a* > /x/T almost everywhere (point 1 of 
Lemma 2.5), then a* can be taken as an admissible signal in (16). 

By homogeneity of the linear system (16) and because i* > 2r, we have that 

V(x(^;0,x*,K,a*)) < 1/2. 

Therefore, for I large enough x(tf, 0, xqj, K, ai) G B\ contradicting (19). ■ 

Let us now turn the general case where the spectrum of A has non-positive real part. The main 
technical difficulties in order to adapt the proof of Theorem 3.1 come from the fact that A may have 
several Jordan blocks of different sizes. 

Theorem 3.2 Let (A,b) G M n (R) x R n be a controllable pair and assume that the eigenvalues of A 
have non-positive real part. Then, for every T,fi with T > jjl > there exists a (T, /j,)- stabilizer for 
(6). 

Proof. Fix a controllable pair (A, b) G M n (R) x R n . Up to a linear change of variable, A and b can 
be written as 

A x A 3 

(n-n')xn' ^2 

where n' G {0, ...,n}, A± G M„/(R) is Hurwitz and all the eigenvalues of A2 G M n _ n /(R) have 
zero real part. From the controllability assumption, we deduce that (^2)^2) is controllable. Setting 
x = (x^ ,x^) T according to the above decomposition, system (1) can be written as 

±1 = A\Xi + ^3X2 + a(t)b\u, (20) 
±2 = A 2 x 2 + a(t)b 2 u. (21) 

If there exists a (T, //)-stabilizer K2 for (21), then 








), h 







K 
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is a (T,(i) stabilizer for (1). It is therefore enough to prove the theorem under the extra hypothesis 
that all eigenvalues of A lie on the imaginary axis. 

Denote the distinct eigenvalues of A by Huj, j G {jo, jo + !>•••> h}, where jo = 1 if cr(A) 
and jo = with too = otherwise. For every j G {0, . . . , h}, let rj be the multiplicity of iuj, with 
the convention that tq = if G" c(-<4). 

Assume that A is decomposed in Jordan blocks. Since (A, b) is controllable, then A has a unique 
(complex) Jordan block associated with each {iuij, — iuij}, jo < j < h. (Otherwise, the rank of 
the matrix (A — iuijld n \ b) would be strictly smaller than n, contradicting the Hautus test for 
controllability.) Therefore, for every j = 1, . . . , h, the Jordan block associated to iujj is ujjA^ + J^, 
where A® = diag(A), . . - , A)) € M 2rj (R) and jf. G M 2rj (R) is defined as 

/ 02x2 Id2 02x2 02x2 \ 

02x2 02x2 Id2 02x2 - - - 02x2 



j: 



c 







2x2 







2x2 



02x2 

Id 2 



02x2 02x2 / 



02x2 
\ 02x2 

that is, in terms of the Kronecker product, = J r . ® Id2- 

All controllable linear control systems associated with a pair (A, b) that have in common the 
eigenvalues of A, counted according to their multiplicity, are state-equivalent, since they can be 
transformed by a linear transformation of coordinates into the same system under companion form. 
We exploit such an equivalence to deduce that, up to a linear transformation of coordinates, (1) can 
be written as 

= Jr xo + ab°u, 
Xj = (ujAW) + J^Xj + aWu, for j = 1, . . . , h, 
where b° and V are respectively the vectors of R r ° and R 2 ^ with all coordinates equal to zero except 
the last one that is equal to one. Here xq G R r ° and Xj G R 2rj for j = 1, . . . , h 

Write the feedback law as u = -K T x = -K^x - Ya=i k I x i with K o G R r ° and Kj G R 2r ^ for 
every 1 < j < h. 

For every u > consider the following change of time-space variables: let 
2/o (*) = D~^ )V x Q {ut), 



(22) 



Vj{t) = (D^^e-^x^t), foil<j<h, 



where D rQtU is defined as in (12) and 



D 



c 



D r ^ u ®Id 2 G M 2rj (R). 
10 



In accordance with 



v(D° tV )- 1 J?.D° tV = J?., D° tV V=V 



we end up with the following linear time-varying system 



< 



ilo = Jr yo - «, 



V 



v 



yi), for j = 1, 



(23) 



where = vD ro> „K , K jyU = vD° Kj and = e -^^ u V for j = 1, . . . , /i. Given v > 0, 

(7) admits a (T, //)-stabilizer if and only if (23) admits a (T/v, /z/i>)-stabilizer. 
For each Z = 1, . . . , h, assume that Kj is of the form (0, k\, . . . , 0, ), that is, 



For uniformity of notations, we also write /Co = Kq ■ 

Let (a u ) u> Q be a family of signals satisfying a u € Q(T/u, fJ,/u) for every ^ > 0. Consider a 
sequence (f n )eN gohrg to infinity as n — > oo such that the matrix-valued curve defined as in 

(9), has a weak-* limit as n — > oo in L°°(R>o, M2h+i-j (R))- Denote the weak-* limit by C*. It 
follows form point 2 of Lemma 2.5 that C*(i) is symmetric and 



for almost every t > 0, for some positive scalar £ only depending on T, /i and cr(.A). 

Define the 2x2 time-dependent matrices Cji, 1 < j,l < h, the 1x2 time-dependent matrices 
Coj, 1 < j < h, and the scalar time-dependent signal Coo by the relation 



Consider, for every n G N, system (23) with v = v n and ^ = X. All coefficients of the sequence 
of systems obtained in this way are weakly-* convergent as n goes to infinity. The limit system is 



We consider (24) as a switched system depending on K in which the admissible switching laws are all 
the time-varying matrix-valued coefficients Cji obtained from the limit procedure described above. 

In the sequel, we only treat the case where is not an eigenvalue of A. The general case presents 
no extra mathematical difficulties and can be treated similarly. Then system (24) takes the form 



K? = ld ®(0,1), JCi = (k[,...,k l ri ). 



C* — \Cjl)j <j,l<h- 




(24) 



ft 




for j = 1,... 



h. 



(25) 



i=i 
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We also assume that the multiplicities r± , . . . , of the eigenvalues of A form a non- increasing se- 
quence. 

Let us impose a further restriction on the structure of the feedback K. Assume that there exist 
k±, . . . ,k ri 6 R, each of them different from zero, such that 

k^ = k ri +i-£, for 1 < I < h and 1 < £ < r;. 

We find it useful to provide an equivalent representation of system (25) in a higher dimensional 
vector space, introducing some redundant variables. In order to do so, for I £ {1, . . . ,ri}, associate 
to y = (yi, . . . , yh) the 2h- vector 

( (/c 1 ^id 2 )(j^y- i 2/1 N 

Yi= \ 

K (/C, i ®Id 2 )(J r ( ;)'- 1 y h / 

Notice that the last 2h — 2m; coordinates of Y\ are equal to zero, where m; denotes the number of 
Jordan blocks of A of size not smaller than I, that is, 

m = #{j | 1 < j < h, Tj > I}. 

For I G {1, . . . ,ri}, let p\ be the orthogonal projection of R 2/l onto R 2m-i x {02/ l _2m i } ) i-e., 

PZ = diag(Id 2m ,, 0(2/ i -2m i )x(2h-2m i ))- 

By construction we have p\ = Id2 ri and p(Yj = Yj for 1 < Z < j < r±. 

Notice that the map (yi, ■ ■ ■ ,yh) l— ► (^l? • • • >^ri) is a bijection between R™ and the subspace 
Et, mr , ofR 2hri defined by 



^....^ = • • • ,Y ri ) | G R 2h and p^ = Y z for Z = 1, . . . ,n}. 

Indeed, the matrix corresponding to the transformation is upper triangular, with the fcj's as elements 
of the diagonal, if one considers the following choice of coordinates on -Emi,...,m ri : take the first two 
coordinates of the first copy of H 2h , then the first two of its second copy and so on until the rf 1 
copy; then take the third and fourth coordinates of the first copy of H 2h and repeat the procedure 
until its rf 1 copy; and so on, until the last two coordinates of the r^ h copy of H 2h . 

If y is a solution of system (25), then Y = (Y\, . . . ,Y ri ) is a trajectory in E^ ^^ satisfying 
the system of equations 

Y^Y^-kPiC^, for Z = l,...,n, (26) 

where, by convention, Y Tl+ \ = 2 ^. 

We prove in the following proposition that there exist k±, . . . ,k n ^ such that system (26), 
restricted to Emi,...,mr > * s ex P° nen tially stable uniformly with respect to all time-dependent mea- 
surable symmetric matrices C* satisfying ^Id 2 ^ < C±{t) < ld2h almost everywhere. 
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Proposition 3.3 For every h, r± G N, for every non-increasing sequence of non-negative numbers 
mi, . . . , m ri such that mi < h and for every £ > 0, there exist A, ki, . . . , k n > and a symmetric 
positive definite 2hri x 2/iri matrix S such that, for every C± G L°°(R>o, M2^(R)) ; if C*(t) is 
symmetric and satisfies £Id2h < C±{t) < ld2h almost everywhere, then any solution Y : R>o — ► 



E: 



of (26) satisfies for almost every t > the inequality 



d 



^-(Y(t) T SY(t))<-X\\Y(t)f. 

Proof. The proof is similar to that of [11, Lemma 4.0] and goes by induction on r\. 

We start the argument for r\ = 1, with h G N, < mi < h and £ > arbitrary. In that case the 
system reduces to 

Yi = -hviC.Yi, 

with Yi G E^ = R 2mi x {02h_2mi}- The conclusion follows by taking ki = 1 and S = ld.2h- 

Let r\ be a positive integer. Assume that the proposition holds true for every positive integer 
j < T\ and for every h G N, < mi < ■ ■ ■ < m ri < /i and £ > 0. Consider system (26) where £ runs 
between 1 and n + 1. Set Y = (Y 2 T , . . . , Y r ^ +1 ) T . Note that if (Yf, . . . , Y r ^ +1 ) T G £t lr .,m ri+1 . then 
Y G ^ 2i ..., mri+1 . The dynamics of (Yi, Y) are given by 



Yi = -kiCYi + UiY, 
Y = -KC*Yi + JY, 



where 



III — (W2/1, 2 / lX 2/i(ri-l))) 



A" 



&2j>2 



J = J ri ® H 2 h. 



Define the linear change of variables (Z\,Z) given by 



Z X =Y U Z = Y + nY u 



where 



n 



( V2P2 



\ r) ri +ip ri +i J 
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and the 77/ 's are scalar constants to be chosen later. Note that Z belongs to E } L n ™ , , if Y does. 
The dynamics of (Z\,Z) is given by 

z x = (-he, + u 1 n)z 1 + ilz, 
z = -((k + + ( J + raii)n)Zi + (J + oiljz. 

Let us apply the induction hypothesis to the system 

Z = {J + nU l )Z, (28) 

which is well defined on -E , m 2 ,...,m ri+ i an d nas the same structure as system (26). (Here C* = Id2h 
and therefore one can take as £ any positive constant smaller than one.) We deduce the existence of 
A > 0, r)i < 0, 2 < I < n + 1 and a symmetric positive definite matrix S such that V"(t) < — A||Z(/j)|| 2 
where V(t) = Z(t) T SZ(t) and Z(t) is any trajectory of (28) in -Em 2 ,...,m ri+ i- Therefore, 

Since O is fixed, for every &i > there exists a unique K(k~i) such that K(fei) + k\^l = 02nhx2h- 
Assume that K = K(k\) and notice that the corresponding &2, . . . , k ri +i are positive. 

Choose S' = (l/2)diag(Id2?i, S) and define the corresponding Lyapunov function W(Z\,Z) = 
\\Zi\\ 2 /2 + Z T SZ/2. If (Zi, Z) is a trajectory of (27), then 

^w(z 1 ,z) = -zf((k 1 c,-u 1 n)z 1 -u 1 z)-z T s((j + nu 1 )nz 1 -(j + nu 1 )z) 

< zfi-hc* + ilo)Zi - A||z|| 2 + (||ni|| + \\s(j + nni)n||)||Zi||||z|| 

< (-fci£ + <5i)||Zi|| 2 - A||Z|| 2 + 5 2 ||Zi||||Z||, 
where the constants 61, 82 > do not depend on k\. Since 

||Zi||||Z|| < e 2 ||Zi|| 2 + " 



s 2 



for every e > 0, then 

^W(Z X ,Z) < f-k^ + Sx + ^j \\Z 1 f + (-\ + e 2 8 2 )\\Z\\ 2 . 
Choosing s 2 small enough in order to have —A + e 2 5 2 < —A/2 and k\ large enough, we have 

±W(Z 1 ,Z)<-^(\\Z 1 f + \\Zf). 

The proof is concluded, since (Zi,Z) and (Y±,Y) are equivalent systems of coordinates on the 

s P ace E m 1 ,...,m rl+ i- ■ 
The proof of Theorem 3.2 is completed by applying the same contradiction argument as in the 
proof of Theorem 3.1. I 
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4 Maximal rates of exponential convergence and divergence 

Let (A, b) £ M n (R) x R n be a controllable pair, K belong to R n and T, p be positive constants such 
that T > jjL. For a € Q(T,p) let \ + {a,K) and \~(a,K) be, respectively, the maximal and minimal 
Lyapunov exponents associated with x = (A — abK T )x, i.e., 

x + / v\ v log(||x(t;0,x , K,a)\\) . . log(\\x(t;0,x ,K,a)\\) 

X r (a,K) = sup iimsup , A (a,K) = ml limml . 

||ao||=l *-»+oo t IM|=1 *-*+°° t 

The ra£e of convergence (respectively, the rate of divergence) associated with the family of systems 
x = (A — abK T )x, a G Q(T,p), is defined as 

vc(A, b, T, p, K) = — sup X + (a,K) (respectively, id(A, b, T, p,, K) = inf X~(a,K)). (29) 



Notice that 



and 



rc(A b, T, u, K) < min mm{-dt(a(A - abK 1 ))}, (30) 
ae[/VT,l] 

rd(A 6, T, u, if) < min min{B((7(A - abK T ))\. 

ae[M/T,l] 

Moreover, since a linear change of coordinates a/ = Px does not affect Lyapunov exponents, then 
rc(A, b, T, p, K) = ic{PAP~ l ,Pb, T, p., {P~ l ) T K) 7 (31) 

and 

rd(A b, T, p., K) = i&{PAP~ l ,Pb, T, p, {P^ fK). (32) 
Define the maximal rate of convergence associated with the PE system x = Ax+abu, a £ G(T, p), 



as 



RC(A,T,p) = sup ic(A,b,T,fi,K), (33) 
A'eR n 

and similarly, the maximal rate of divergence as 

KD(A,T,p,) = sup rd(A,b,T,p,,K). (34) 
KeR n 

Notice that neither RC(A, T, p) nor KD(A,T, p) depend on b, as it follows from (31) and (32). 

Remark 4.1 Let us collect some properties of RC and RD that follow directly from their definition. 
First of all, one has 

RC(yl + AId„, T, p) = RC(A, T, p) - A, RD(^4 + AId n , T, p) = RD(A, T, p) + A. (35) 
Then, by time-rescaling, 

RC(A T, P T) = RC(A/T, 1, p), RD(A, T, pT) = KD(A/T, 1, p). (36) 
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Notice moreover that, thanks to (13), both RC(J n ,T, pT) and RD(J n ,T, pT) only depend on p and 
thus are equal to RC( J n , 1, p) and RD( J n , 1, p), respectively. Finally, because of point 2 in Lemma 2.4, 
RC and RD are monotone with respect to their third argument. 

Remark 4.2 Given a controllable pair (A,b) and a class Q(T,p) of PE signals, whether or not RC 
and RD are both infinite can be understood as whether or not a pole-shifting type property holds 
true for the PE control system x = Ax + abu, a G G(T, p). 

The study of the pole-shifting type property for two-dimensional PE systems actually reduces to 
that of their maximal rates of convergence as a consequence of the following property. 

Proposition 4.3 Consider the two-dimensional PE systems x = Ax + abu, a G Q(T,p), with (A, b) 
controllable. Then RC(^4, T, p) = +oo if and only ifRD(A,T,p) = +oo. 

Proof. According to (31), (32) and (35), it is enough to prove the result for (A,b) in companion 
form and with Tr(^4) = 0. Let then 











\ -i 




[a 1 







A = 
with a E R. 

Assume that KC(A,T, p) = +oo. By definition, for every C > there exists K G R 2 such that 
rc(>l, 6, T, p, k) > C. Therefore, by definition of rc, 

limsup log(||x(t;0,x ,^a)||) < _^ ^ g g^^y^ = L (38) 

t^+oo t 

Moreover, due to (30), for C large enough we can assume that k\,k2 and k\lki are large positive 
numbers. 

Let K_ = (ki,—k2). We claim that if C is large enough then RD(A, b, T, p, K_) > C. Assume 
by contradiction that there exists a G Q(T,p) such that A~(a,i^_) < C. Then there exists x G R 2 
of norm one and an increasing sequence (£ n )neN of positive times going to infinity such that 

log(\\x(t n ;Q,x,K-,a) 



Notice that for every t G [0,£ n ], 



< C, V G N. 



x(t;0,x,K_,a(-)) = diag(l, -l)x(t n - t;0,x n ,K,a(t n - •)), 



where x n = diag(l, — l)x(t n ; 0, x, a) 
Therefore, by homogeneity, 



log 



x (i n ;0, T^u,K,a(t n 



lo g(lkn||) = log(\\x(t n ;0,x,K_,a)\\) > _^ 
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This would contradict (38) if, for some positive integer n, x n /\\x n \\ = x and the signal obtained 

by repeating <5|[o,t„) by periodicity over R>o belonged to Q(T,fj,). Indeed, in such a case, 

log(||a(fct w ;0,x,ir,a(-))||) > _ c 
kt n 

for every k > 1, where a G Q(T,fj,) denotes the signal obtained by repeating a\^ ^(t n — •) by 
periodicity over R>o- 

In order to recover the periodic case, we are going to extend a backwards in time over an 
interval [— 2[i — T n ,0) as follows. First set A± = A — bK-£. We take a = 1 on the intervals 
[— fi, 0) and [—2fi — T n ,—fi — r — n) and we extend a on [— fj, — r n ,—fi) in such a way that the 
trajectory corresponding to a|[_ /J „ Tn _ M ) and to the gain connects the half-line R>ox„ to x~ , 
where x+ = exp(/^4j~)diag(l, — l)x n and x~ = exp(— fiA^)x. We show below that this can be done 
fulfilling the PE condition and with r n upper bounded by a constant independent of n. Hence, the 
signal obtained extending a[_ 2At _ Tni t n ] by periodicity belongs to C?(T, /i) and we have 

x(t n + 2fi + T n ;0,x n ,K, a(t n + 2fi + T n - •)) G R>o£ n 
log ^ x ^t n + 2fjL + T n ;0,-^j,K,a(t n + 2fi + T n - •)) ^ = log (||x||) - log(||x(t n ; 0, x, a)\\), 

where x = x(r n + 2/i; 0, diag(l, — l)x, K, a|[_ 2M _ rni o] ( — '))■ Note that log(||x||) can be lower bounded 
independently of n, because of the uniform boundedness of r n . Therefore, 

log(\\x{t n + 2n + T n ;0,x n ,K,a{t n + 2^ + T n - -))||) log (p||) Ct n 



t n + 2fi + T n t n + 2fl + T n t n + 2fl + T n 

is larger than — C for n large enough and we can conclude as in (40). 

We are left to prove that the control system on the unit circle whose admissible velocities are 
the projections of the linear vector fields x ^ (A — ^bK^)x, £ G [0,1], is completely controllable 
in finite time by controls £ = £(t) satisfying the PE condition. Notice that the equilibria of the 
projection of a linear vector field x \— ► A'x on the unit circle are given by the eigenvalues of A'. All 
other trajectories are heteroclinic connections between the equilibria, unless the eigenvalues of A' are 
non-real, in which case the phase portrait is given by a single periodic trajectory. 

Denote by 9 a point on the unit circle, identified with R/2-7rZ. Then, the above mentioned control 
system on the unit circle can be written 

9 = acos 2 (e)-sm 2 (e)+^cos(e){k2sm(e)-k lC os{9)), £g[0,1]. (41) 

We prove the controllability of (41) by exhibiting a trajectory 9 of (41) corresponding to a PE 
control £, starting at some 9q G R/2-7rZ, making a complete turn and going back in finite time to #o- 

The PE condition will be verified by checking that the control £ = is applied for a total time 
that is smaller than T — jjL. Define the angle 9k G (0, vr/2) by 

ko 

tan. (0 K ) = 2^. 

h 
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Notice that the eigenvectors of A^/ are proportional to the vectors (2, k% =t ^/k% — 4(/si — a)). There- 
fore, assuming that k\ is larger than a, the angle between any real eigenvector of and the vertical 
axis is smaller than Ok- 

Take Oq = tt/2 and apply £ = until reaches tt/2 — Ok- Since A^/^i is small and Ok is of the 
same order as fo/fci, then we can assume that acos 2 (#) — sin 2 (#) < —1/2 for 9 £ [tt/2 — Ok, 7T /2]. 
Therefore, the time needed to go from tt/2 to tt/2 — Ok can be assumed to be smaller than (T — p)/2. 
When the trajectory reaches tt/2 — Ok, switch to £ = 1 and apply it until reaches (in finite time) 
—tt/2. This is possible since either the eigenvectors of are non-real or they are contained in the 
cone 

{(rcos0,rsin0) | r > 0, G (tt/2 - K + rmr, tt/2 + mi), m G Z}. 

In both cases the dynamics of (41) with £ = 1 describe a non-singular clockwise rotation on the arc 
of the unit circle corresponding to [tt/2, tt/2 — Ok]- The trajectory is completed, by homogeneity, 
taking £ = until reaches —tt/2 — Ok and finally £ = 1 until reaches — 3tt/2 = tt/2 (mod 2tt). As 
required, the sum of the lengths of the intervals on which £ = does not exceed T — p. 

This concludes the proof that RC(^4, T, p) = +oo implies RD(A, T, p) = +oo. The converse can 
be proven by a perfectly analogous argument. ■ 

4.1 Arbitrary rates of convergence and divergence for p large enough 

This section aims at proving that for p large enough a persistently excited system can be either 
stabilized with an arbitrarily large rate of exponential convergence or destabilized with an arbitrarily 
large rate of exponential divergence. This will be done by adapting the classical high-gain technique. 

Proposition 4.4 Let n be a positive integer. There exists p* G (0, 1) such that for every controllable 
pair [A, b) G M„(R) x R", every T > and every pe(p*, 1] one has RC(A, T, pT) = RD(A, T, pT) = 
+oo. 

Proof. Fix T > and let (^4, b) G M n (R) x R™ be a controllable pair in companion form. According 
to (35), it is enough to establish the result with the extra hypothesis that Tr(A) = 0. We therefore 
assume in the sequel that b = (0, . . . , 0, 1) T , A = J n + bK^ and K^b = 0. 

We first prove the stabilization result. Fix K G R n such that J n — bK T is Hurwitz. Let P be 
the unique positive definite n x n matrix that solves the Lyapunov equation 

{J n - bK T ) T P + P{J n - bK T ) = -Id n . 

Define V(x) = x T Px. Then, for every a G L°°(R, [0, 1]) and every solution of x = (J n — abK T )x, 
one has 

jV(x(t)) < -dVixit)) + C 2 (l - a(t))V(x(t)), 

IS 



with Ci, Ci two positive constants only depending on K. Choose p £ (0, 1) and assume that a is a 
(T, Tp)-signal. Then, for every t > 0, 

^(x(t + T)) < F(x(i)) exp(-r(C 1 - C 2 (l - p))). 

Therefore, if p > 1 - (C 1 /2C 2 ) then RC(J n ,T,Tp) > Ci/2 > 0. For every 7 > 0, set K 7 = jD 7 K 
(where, as in the previous section, D 1 = diag(7 n_1 , ... ,7, 1)). Recall that J n and L> 7 satisfy (13). 
Take a solution of x = (A — abK^)x with a £ Q(T,pT). Set z(-) = D^x(-) and notice that for every 
7 > 1 

jV(z(t)) < 7(-Ca + C 2 (l - a(t)) + ^/ 7 2 )^(z(t)), 

where only depends on and P. Then clearly RC(A,T,Tp) > jCi/3 for p > 1 — (C1/2C2) 
and 7 large enough. Thus, RC(j4, T, Tp) = +00 and one can choose p* > 1 — (C1/2C2). 

The destabilization result can be obtained by a similar argument based on the Lyapunov equation 

{J n - bL T ) T Q + Q(J n - bL T ) = Id n , 

verified for some L £ R n and some symmetric positive definite matrix Q. ■ 

4.2 Finite maximal rate of convergence for p small enough 

In this section we restrict our attention to the case n = 2. 

Proposition 4.5 There exists p* £ (0, 1) such that for every controllable pair (A, b) £ M2(R) x R 2 , 
every T > and every p £ (0, p*) one has KC(A,T, pT) < +00. 

Proof. Thanks to Remark 4.1, it suffices to show that there exists p* £ (0, 1) such that, for every 
controllable pair (A,b) £ M 2 (R) x R 2 with Tr(A) = 0, one has RC(A, l,p*) < +00. 
As in (37), take (A, b) in companion form, ie, 

A = J 2 + aH, 6=(0,1) T , 



with a £ R and if = 




For £ [— 7r,7r) set eg = (sin#, cos 6) T and define 7/0 = ( — 1> 0) T - Every gain can be written as 

K e n = 7- D 7 e 6», 

with 7 > and 9 £ [— 7r,7r). 

Moreover, if ^4 — 6i^ T is Hurwitz with X = ^jD^eg then the sum and the product of its two 
eigenvalues are, respectively, 7cos# > and 7 2 sin# — a > 0. In particular, 9 £ (— tt/2, it/2) and 
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7 2 sin# > a. If 9 G (—ir/2, 0] with A — bK T Hurwitz, then \a — sin#7 2 | < \a\ = —a and therefore the 
convergence rate of A — bK T is upper bounded by a constant only depending on a. 

Let Slo = (0, vr/2) x (0,oo). We show in the following the existence of p > and Q = {(#,7) | 
< 9 < vr/2, < 7 < 7(0)} C fi such that 

if (9, 7) G Slo and Kg ~ is a (1, /^-stabilizer of £ = Ax + a&u, then (0, 7) G fi, (42) 

and 

sup mm{-SR(a(A - bKj ))} < +co, (43) 
(e,7)en 

and the conclusion then follows from (30). 

Fix 9 G (0,7r/2). In order to find, for 7 large enough, a G £7(1, p) and Xo G R 2 such that the 
trajectory of 

x = Ax — abKg^x, x(0) = xo, 

is unbounded, we apply the transformation y 7 (-) = Djx(-/j): the problem is now to find, for 7 large 
enough, a G £7(7, p'j) and an unbounded trajectory of 

V = (j2 + ^H^j y - abegy. (44) 

Due to the homogeneity of the system, the latter fact reduces to determine r large enough and 
a G £7(t, 2pr) such that the solution y{- ; 0, yo, eg, a) of (44) satisfies y(r; 0, yo,a) = —£2/0 with £ > 1. 
Indeed, for every 7 > r the extension of a|[o, r ) by periodicity is a (7, p7)-signal (see point 4 in 
Lemma 2.4) and the sequence ||y(mr; 0, yo, a)\\ = £ m goes to infinity as m goes to infinity. 

Set 

M e = J 2 - bej, N afia = J 2 + - bej . 

Consider h > small to be fixed later. We distinguish two cases depending on whether 9 G (0, h) 
or not. 

The case 9 G [h,ir/2). 

We construct a PE signal a as follows: starting at yo take a = 1 until the trajectory y(- ; 0, yo> e e, 01) 
of (44) reaches, at time T\, the switching line sin(#)x + cos(9)y = 0. In order to ensure that the 
switching line is reached in finite time and, moreover, that T\ is lower and upper bounded by two pos- 
itive constants only depending on h (and not on 9 G [h, vr/2)), it suffices to choose 7 > Ti(a,h) > 
with Ti(a, h) only depending on a and h. (Indeed, the bounds hold for all matrices in a neighborhood 
of {Mg I 9 G [h, it/2)} and it suffices to ensure that N a< g^ belongs to such neighborhood.) 

From y(Ti; 0, yo, eg, a) set a = until the first coordinate of ?/(• ; 0,yo, eg, a) takes, at time 
T\ + T2, the value 1. Finally, take a = 1 until the second coordinate of y(- ;0,yo,eg,a) reaches, at 
time T1 + T2 + T3, the value 0. (See Figure 1.) 
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\ sin(9)x + cos(0)y = 
















i » 






-Vo 


Svo 



Figure 1: The trajectory y(- ; 0, yo, eg, a) when G [/i, vr/2) 



Analogously to what happens for T\, the values T2 and T3 admit lower and upper positive bounds 
only depending on h. 

Define r = T\ + T2 + T3 and notice that it admits an upper bound T\ (h) only depending on h. 
Finally, T ^^^ T admits a lower bound p\ only depending on h. The construction of the required 
(r, pir)-signal is achieved and we set 

j(e) = max(T 1 (a,h),T 1 (h)). (45) 



The case 9 G (0,/t). 

Notice that the condition for N a g~ to be Hurwitz is that 7 2 > |a|/sin#. Choose 7 > ^(a, #) = 
Myjof/ sin# with M large (to be fixed later independently of all parameters). In particular, for M 
large enough and ho > small enough (independent of all parameters), for every 9 € (0, /to) and 
every 7 > r2(a,0) the matrix N a< g^ has two real eigenvalues, denoted by {J>+(a, 9,j) > fJ,-(a, #,7) 
and 

-2 < /i_(a,0,7) < -1/2, -2 sin < ^+(0,^,7) < -sin 9/2. (46) 

From now on we assume h € (0, /io). 

Similarly to what has been done above, we construct a PE signal a as follows: starting at 
yo take a = 1 in (44) for a time Ti = pM/\fi + (a, 9, 7)! with /) G (0,1) to be fixed later. Set 
Vi = y(Tr,0,yo,eo,a). 

From yi set a = for a time X2 = M/\fi + (a, 9, 7)! and denote by 7/2 the point y(Ti+T2; 0, j/0j eg, a). 
Finally, take a = 1 until the second coordinate of y(- ; 0, yo> &g,a) assumes, at time T\ + T2 + T3, the 
value 0. (See Figure 2.) 

We next show that there exist p and M independent of 9 and a such that T3 is well defined and 
y(Ti + T 2 + T 3 ; 0, y , e e , a) = -£y with ^ > 1. 
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-3/o -t yo 



Figure 2: The trajectory y(- ; 0, yo, eg, a) when G (0, h) 



A simple computation yields 

1 

yi 



H-(a,9,j) -fi + (a,0,j) \ fl _(a,e,j)fi + (a,9,-f)(e^- ( - a ' e ^ Tl - e ^ a ' 9 '^ Tl ) 



-pM 



-1 



+ 0(9 2 ), 



with ||0(# 2 )|| < CO 2 and C only depending on M and p. (Similarly, in the sequel the symbol O(0) 
stands for a function of 8 upper bounded by CO with C only depending on M and p.) 
In addition, one also gets that the first coordinate of yi is equal to 



e~ M P(M -1)+O(0) 

M p ( M i^2) sinh ( _^n» 

Af p / M afep sin C sin 



cosh 



cos 



if a = 0, 

S me ^+0^ if a <o. 



^ + (a,0,7) 

Using (46) one deduces that the first coordinate of 2/2 is larger than 

e -^P(M/2sinh(l/2) - cosh(2)) + O(0) if a > 0, 
e -Afp(M/2sin(l/2) - cos(2)) + O(0) if a < 0. 

Then in all three cases the first coordinate of j/2 becomes larger than 

e- M P(MC -C 1 + O(0)), 

and one also gets that the second coordinate of yi can always be lower bounded by 

sin0e- M ^(<7i - C /M + O(0)), 

with Co > and C\ > independent of all the parameters. 
Fix M large and p G (0, 1) such that 



-Mp 



(MC - Ci) > 2, e - A ^(Ci - C /M) > Ci/2. 



Finally, by eventually reducing /i in order to make each O(0) uniformly small, one can ensure 
that the first coordinate of y2 remains larger than 1 and that its second coordinate is positive. 



22 



Similar computations to the ones provided above show that it is possible to further ensure that 
T 3 < 2Ti. 

Define t = T 1 +T 2 + T Z . Then M/(2 sin 9) < r < 8M/ sin 9 = T 2 (9). Choose now 

7(0) =Af(8 + v / ^)/sin6>>max(T2(fl),r 2 (a,60). (47) 
By construction, a G G(t,p~t). To conclude the proof it is enough to check condition (43) on 

= {(0, 7 ) | < < /t, 0<7< 7(0)}. 



For (0, 7) G fi* define 



= A - bKle 



1 

a — 7 2 sin 9 —7 cos ( 



Then 



< det(A^ b ) < C |Tr(A^ b )| + \a\, 
with C = 2M(8 + ^/faj), implying (43). ■ 
The following corollary is a direct consequence of Remark 4.1 and Proposition 4.5. 

Corollary 4.6 Take p* as in the statement of Proposition 4-5- For every controllable pair (A,b) G 
M 2 (R) x R 2 , every T > and every p < p*, if A > is large enough, then (A + AId2,6) is not 
(T, pT)-stabilizable. Moreover, if < p < p* and A > RC(J2,l,p), then (J 2 + AId2,6o) * s n °t 
(T, pT)-stabilizable for every T > 0. 

The above corollary establishes the existence of non-stabilizable PE systems if the ratio p = p/T > 
is small enough and regardless of T. This is rather intriguing when one recalls, on the one hand, 
that any weak-* limit point a* of a sequence (a n ), with a n G Q(T n , pT n ) and lim n _ >+oo r n = 0, 
takes values in [p, 1] (see point 1 of Lemma 2.5) and, on the other hand, that the switched system 
x = J 2 x + a+fybou, u+(t) G [p, 1], can be uniformly stabilized with an arbitrary rate of convergence 
by taking the feedback law u-y = —jD.yKx, where 7 > is arbitrarily large and K is provided by 
[11, Lemma 4.0]. 

Remark 4.7 One possible interpretation of Proposition 4.5 goes as follows. Consider the desta- 
bilizing signals built in the argument of the proposition back in the original time-scale, i.e., as 
(1, p) -signals. These signals take only the values 0,1 over time intervals of length proportional to 
I/7. Therefore, the fundamental solution associated to x = (A — ab^K^^x is a power of the product 
AiA 2 Az, where A ± = exp(ri(A - b K ltd )/j), A 2 = exp(T 2 A/7) and A 3 = exp(T 3 (A - b Q K lfi )/^). 
The stabilizing effect of A — b^K^ g is countered by the overshoot phenomenon occurring when the 
exponential of A — boK^ g is taken only over small intervals of time. If 7 is large enough, such 
overshoot eventually destabilizes x = (A — aboK^^)x. 
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4.3 Further discussion on the maximal rate of convergence 

Let (A, b) £ M(n, R) x R n be a controllable pair. Define 

p{A, T) = inf{p £ (0, 1] | RC(A T, Tp) = +oo}. (48) 

Notice that p(A,T) is equal to p(A/T, 1) and does not depend on Tr(^4) (see Remark 4.1). 

Proposition 4.4 implies that p(A, T) < p* for some p* £ (0, 1) only depending on n. In the case 
n = 2, moreover Proposition 4.5 establishes a uniform lower bound p(A,T) > p* > 0. 

The following lemma collects some further properties of the function T \— > p(A,T). 

Lemma 4.8 Let (A,b) £ M n (R) x R n be a controllable pair. Then (i) T ^ p(A,T) is locally 
Lipschitz on (0,+oo); (ii) there exist limr^+oo p(A, T) = sup T>0 p(A, T) and lim T p(A,T) = 
mi T>0 p(A,T). 

Proof. In order to prove (i), notice that point 3 in Lemma 2.4 implies that if RC(A,T, pT) < +oo 
then for every rj S (0, pT), 

RC (a,T +v ,^-(T + r } )\ < +oc, (49) 

RC (a, T-t), ^Y^j-(T - rj)) < +oo. (50) 

From (49) we deduce that for every rj £ (0, p(A,T)T), 

p(A,T + v )> p -^p^, (51) 
T + rj 

and thus 

p(A,T)-p(A,T + rj)<rj/T. 
Similarly, (50) implies that, for every rj € (0, p(A, T)T), 

1 — rj 

Therefore, one has 

MA,T)> * AT + lf + "»-" (52) 

for every r\ satisfying < rj < p(A, T + rj)(T + rj)) and in particular for every rj £ (0, p(A, T)T) (see 
(51)). We obtain from (52) that p(A,T + rj) — p(A,T) < n/T and we conclude that 

\p(A,T + V )-p(A,T)\<^ 

for every rj £ (0, p(A, T)T). 

As for point (ii), it suffices to deduce from point 5 in Lemma 2.4 that if < p' < p < 1 then 
there exists M > such that whenever RC(A, T, pT) = +oo one has RC(^4, 7, p'j) = +00 for every 
7 > such that j/T > M. ■ 



24 



Remark 4.9 In the case A = J n equality (15) implies that the function T *— > p(J n ,T) is constant. 
When n = 2 its constant value is positive, due to Proposition 4.5. 

5 Open problems 

We conclude the paper by providing some questions that arose from our investigation of single-input 
persistently excited linear systems. 

Open problem 1 Does Proposition 4.3 still hold true in dimension bigger than two? Notice that 
the proof provided here essentially relies on the controllability of (41) in finite time. 

Open problem 2 Consider the constant p* defined as the upper lower bound for all the p*'s satis- 
fying the statement of Proposition 4.4 (n fixed). What can be said on the dependance of p* on n as 

n — ► oo? 

Open problem 3 We conjecture that Proposition 4.5 holds true in dimension n > 2. Note however 
that the proof given in the 2D case cannot be easily extended to the case in which n > 2. Indeed, 
our strategy is based on a complete parameterization of the candidate feedbacks for stabilization 
and on the explicit construction of a destabilizing signal a for every value of the parameter 8, which 
takes values in the one-dimensional sphere. In the general case, the parameter would belong to an 
(n — l)-dimensional manifold and an explicit construction, if possible, would be more intricate. 

Open problem 4 It is a challenging question to determine whether the function T i— > p(A, T) 
(defined in (48)) is constant for a general matrix A. If this is true, one may wonder whether its 
constant value depends on A. Otherwise, a natural question would be to understand the dependence 
of lim^_^o+ p(A, T) and limr^+oo p(A, T) on the matrix A. 

Open problem 5 Proposition 4.5 states that, for n = 2 and p/T small, the PE control system x = 
Ax + abu, a 6 Q(T,p), does not have the pole-shifting property (see Remark 4.2). It makes therefore 
sense to investigate additional conditions to impose on the PE signals (periodicity, positive dwell- 
time, uniform bounds on the derivative of the PE signal, etc) so that the pole-shifting property holds 
true for these restricted classes of PE signals, regardless of the ratio p/T. First of all, the subclass of 
periodic PE signals must be excluded, since the destabilizing inputs constructed in Proposition 4.5 
are periodic. It is also clear that, for the subclass of G(T, p) given by all signals with a positive dwell 
time td > 0, one gets arbitrary rate of convergence (or divergence) with a linear constant feedback, 
for every choice of T, p, td- Here follows our conjecture. 

Given T, M > and p G (0,1], let V(T, p, M) be the subset of Q(T,pT) whose signals are 
globally Lipschitz over [0, +oo) with Lipschitz constant bounded by M. Then, given a controllable 
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pair (A, b), we conjecture that it is possible to stabilize (respectively, destabilize) by a linear feedback 
the system x = Ax + abu, a G T>(T, p, M), with an arbitrarily large rate of convergence (respectively, 
divergence), i.e., we conjecture that for every C > there exist two gains K\ and K2 such that for 
every a £ T>(T, p, M) the maximal Lyapunov exponent of x = (A — abKj)x is smaller than — C and 
the the minimal Lyapunov exponent of x = (A — abK^x is larger than C. 
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